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Tropical cyclone recurvature is one of the most 
difficult events to forecast for even the most experienced 
forecaster. The ability to forecast recurvature correctly 
could annually save both the military and civillian 
communities millions cf dollars, “thousands of man hours and 
the lives of many victims. There are £2w reliable rules 
that the Bmapıcalsicyerone Frerecaster can follows eros 
forecasting  recurvature cr non-recurvaturs in either the 
near and the extended time frames. Riehl and Shafer (1944) 
were amorg the first to summarize guidelines for forecasting 
Beeurvature cl Atlantic tropical cyclones. Rishl and Shafer 
defined  recurvature as "any major alteration of a storm 
meaek" which 2s quite non-descri pt. Riehl and Shafer listed 
Bew-real factors affecting tropical cyclones motion, including 


Bae] effects of the upper-level winds and the interaction 


Ben <roughs in the westerly flow `o the west of the 
meopical  cycione. Tacy also GeSGEL DSi Vear2ous Synoptic 


Eta tions to assist the forecaster in predicting the storm 


Eh. The Jcint Typhcon Warning Cente 
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Ileana. "et auchopteal Cyclone Erom an initial p 
toward the west and northwest to the north then northea 


which is also quite non-descript. 

More recent studies by George and Gray (1977), Char, et 
al (1980) and Shaffer (1982) have shown that the tropical 
cyclone motion is largely dependent upon the surrounding 
steering flow represented by height/temperature fieids. 

The purpose of this study is to use a simple 
representation of the synoptic foccing by Epes i Gee. 
Szsehogonel Functions (EOF's) plus persistence-related 
variables tc predict ë tropicil cyclone reenrwasure Or 
Won-recurvature. Tae impe incréðents of 36, 54 and 72 h 
for forecasting recurvaturs were selected fer several 
reascns. Huc aminimumes:he perazod of 36 h was thought 
to be of greatest value and to provide the nost confident 
forecast for the shore establishment, air wing squadrons and 
ships presentiy in port. Secondly, The upper-air soundings 
mee transmitted to Fleet Numerical ‘Oceanography Center 
(NOC) by approximately 3h after 00/12 GMT. These 
soundings are used tc generate the 12-n2urly forecasts of 
mee 500 mb hei-gh=s (in this case). These generated fieids 


Mould be available tc Jcint Typhoon Warning Center (JTWC) 


11 





a 


Guam DNI Oc cxuWwarnsnd at approximately 6 h 
after 00/12 GMT. The computer program to decompose the 
D-value fields (deviations from the standard atonsphere 
gzopotential) into the eigenvectors (EOF coefficients) and 
compute the persistence statistics, would take at most a few 
minutes, and could actually be done on a hand hela 
salentator (Shaffer, 1982). This would make the warning 
time increments equal to 30, 48 and 66 h i£ the forecast 
technique makes use Or the analyzed and forecast 
@=Opotential height fields. This time delay problem could 
be partialiy eliminated by using the previous 12-h forecast 


m elüs to define steering, which would permit the technique 


ct 


to be availabie at the time JTWC prepare Heir forecast es 


Ul 


00/12 GUT. It is believed that 1ittla effect would be 
apparent for the shorter forecast periods, although there 
would likely be a larger impact on the extended time perio 


forecasts. 
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of Brown 
Borcing to 
North Pacific Ocean by 
selection 
1981; 

Exactly 


the fields and 


Eure 1. G 


(1981) and an EOF representation of 


Predscr Broca cyclone motion in 


Shaffər 


of storm criłŁłeria (see dascripticrs 
Shaffer, 1982; and Shaffer and ELsberry, 
the same in all three studies. A short 
eriter oa. fOJILOW. ; 
+ ta Tot tot + + + + + + + + + 225 
+ ++ e+ eee eee e+ + + + + 20 
2 1 L 
ec ae ee + S A 
ee er re ro 1 
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5 13 D 
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40 30 20 10 0 -0 -20 -30 
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uid Arrangement Usei in Statistical 
e an lodi cal Track Prediction. 
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II. EOF ANALYSIS OF SYNOPTIC FORCING APPROACH 
TO” RECURV ATURE 
This research is a sequel to a synoptic map type scheme 


the synopti 


Summary of 





The 500 mb D-values are extracted on a movable 120-point 
mead (fig. 1) Over 35° in latitide and 709 in longitude with 
a 50 nerement. The tropical cyclone is always centered at 
grid point 70 which is 10° from the southernmost boundary 
and 40° from the western boundary. This choice is made to 
include botn a OB Ton Of the nid-latituds flow 
characteristics and asmaller section of the more uniforn 
tropical height fields to the  scuth of the storm. The 
movabie feature is incorporated to minimize the seasonal 
Variation that would occur if the grid was stationary. The 
earth's curvature (Ccnvergence of the maridians) is not 
includei over the grid but could be easily addad for 
cperational use. The D-values on the standazd 63 x 63 grid 
from che FNOC archives are deviations from the standard 
atmosphere geopotential and are fitted to the 120 point grid 
by a Bessel linear interpolation scheme. 

nequezements for -Topical cyclones Ło be included ina 


this study are: 


U} 


Ton cal Cyclone must 
sity (>17 m/sec) an 


Ie yclecns must be locat=d between 10°N ard 250N (This 
Ss done to prevert the grid from Seer eng Le souütheri 
hemisphere and to avoid "possible problems eee EC 
the convergence cr the meridians at higher Latitudes); 
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2 for the scheme to be 


e The storms must exist at 0)/ GMT 
ng FNOC data analyses; 


1 
able to use the corresponding 


e Sach storm observation inclujed must have 35 hours 
5etween a prior observation of the same storm *o insure 
pseudo-independence withia the sample. This ver 
noue crar restrictlion ?nusures that che statistica 
results are unbiased. 


Ore approach to assess spatial variations in geophysical 
fields is through the use of Empirical Orthogonal Function 
(EOF) analysis (also known as Principal Component analysis). 
EOF analysis may be applied to any multi-dimensional data 
set to separate Signal from noise, and identify significant 
Spatial correlation patterns {Hotalling (1933); Kutzbach 


BER: Morrison (1967); Barnett (1978); Richman (1980); 


(D 
(= 
un 
D 
~ 
cr 
E 
(D 
Q 
f» 
(1 
[D 


Shaffer (1982) and others). ror geophysical 

are ovolaced on the moveable grid with sach grid point 
representing a dimension in space. The 12t2 are reoresented 
Bee metrix, wath M the number cf grid peints, and N 
the number oí data cases. Each colunn Ofthħhisi matri: 


represents a single realization of the observed geophysical 


Bela. A Maa 2S them tormed by multiplyiag the 


tJ 


Cata matrix by its transpose. The form of R depends upo! 


Fu 


the ca nthe Initial data matrix, but is gensrally, 


Aer a COVariance ot cocrelation matrix of <the zd 


“Q 
13 


pornts. ir cne wishes to give equal weighting to each arid 


15 





A 





Beem, tne correlation matrix is formed. ' If, on the other 
hand, the desire is to give greatest weighting to those 
dimensions (grid points) with the largest univariate 
variation, the covariance matrix is formei. 

The choice of correlation versus covariance matrix is 
Ber ctical for meteorological fields whic cover a large 
spatial domain. For example, in analyzing geopotential 
fields, mid-latitude qid o o:nts have much greater 
Dur2ability tnan trofical grid points. Inr "he am nis to 
focus on mid-iatitude features, the covariance matrix forn 
of R is preferable. ite Onmrne Other hand, zhe wish iS ze 
determine variability patterns over the entire domain, the 
See-elation natrix form cÉ R is favored. Ries correlation 


matrix approach mutes the EOF response in higher variability 


I 


regions (mid-latitudes), while amplifying the response in 


tN 


(p 
ct 
D 
in 
ct 
(D 
ps 


BR gcons cf lower variability. Both approaches w 


(D 
p 


this study. With the covariance matrix approach, thes 


Pome sxpralzned varietlen in the D-value fields for the 


first 10 eigenvectors wes more than 94%, which was 9% 
Heater than the correlation eigenvectors. Regression 


analyses from the covariance method were shown 0 be 


Brot to zhese from the correlaticn msthod. TRUS, FR 


(D 


ear tr elation matrix metho 


¡Y 
pa 
Ui 
f 
Ui 
(U 
As 
pa 
ej 
ct 
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te 
Ul 
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(2 
fu 
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Micon fOr ng He ma Tis R, minimization of the total 
variability leads Stectivaet Oo the sigenvalue problen 
(Morrison, 1967). The eigsnvector associated with the 
largest eigenvalue is the single vector, out of all possible 
vectors in space, that minimizes the summed variation 
(squared) of all of the observation vectors in M space. 


Furthermore, this eigenmods explains exactly 


of the total variation in the observations, where A is the 
eigenvalue of the ith mode. Since rhe. trace of a 
e elation matrix is identically equal *o the order of the 
mc ix, the percent of total variation explained by any 


eigenvector is simply 





The elements cf each eigenvector may be loosely interpreted 
@emene MULci Die correlation of the correspondirg grid voint 
ER ali other grid points, Simultaneously. The exact 
meaning of the eigenvector elements  (loalings) is given by 


MOr rison (1967). 


17 





The aim of this project is to investigate spatial 
variation structures cf analyz2d 500 mb D-value fields for 
the western North Pacific Ocean. Therefore, we choose the 
correlation matrix approach. fie “yet iave step in data 


processing is to point-normaliz2 ths data. Specifically, we 


use 
N 
aas een Di, 5) 
N j:1 
pg) = DEI - Dba Ete) 
N 1 
b> (05:05 3) e Doar) je | 2 
| j: 1 N -1 


Where j is the map (case) index and i is the grid point 
index. nuc Dx: 3) 2s Simply a no»ruJally distributed 
random variable, With a mean at each grid point over all 
maps of zero and z variance of 1. This allows an 


equally-weighted, geid point interpretation S the 


w 


eigenvectors. Lomm his SUI y, "ens data 


(v 
t 
fv 
pa 

M 
{N 
ya 
Uu 
£ 
n 
N 


Sempactsd using the 500 mb fields only. 


The date corsist of D-values of geopotential, (the 
difference between the observed geopotential and the 
standard gecpotantial, SO that 2 positive D-value 
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Potentials), 


in 


—- 


muecwrtrrst 9 EOF 


2 TE 


corresponds to higher than standard gs 
meters. EOF analysis was performed for 120-dimension space 
E15 grid) in Fig. 1. Examples of 
eigenvectors (correlation case) are shown 


A 


must be stressed that the neaative of the pattern shown in 
these figures is also a solution to the problem. (For 
examples of 700mb and 850mb eigenvectors see Shafter 
(1982) .) 





Eiqenvectors 1-9 
Indicates the Lcc 
Ber dw orte). 


Pegure 2. 


The number of eigenmodes 


determined by compar ing 
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ie Ore. nee Xk 

FOoprcai Cyclone, 

in each case is 
of the analyzed 





geopotential fields with the aigenvalues of a Monte Carlo 
Simulation cf random fields as used by Preisendorfer and 
Barnett (1977). The eigenmodes retained are those that are 
larger than the random eigenvalues plus 3 standard 
deviations. By using this approach we are 99% confident 
that the only eignvsctors ratained arz those containing 
actual signal, and that we are discarding noise only. Using 
this selecticn procedure, the actual nuaber of significant 
eigenvectors (EOF's) retained 183 limited to 10. This number 
of EOFs explains approximately 85% of tha variation ir the. 
perds. Table I outlines tne increase of explained variance 
as the number of significant eigenvectors retained is 
Mijeseesed. rig. 3 is tae graphical repras=ntation of Table 
I. 

The total sampl2 of the 504 storms during the 1967 to 
1976 time period ¡is divided into a dependent sample cf 454 
cases and a 50-case independent sample Zar use as the test 
sanple. Demi aes dent: cal Sampis of 50 randomly 
Selected independent storms that Shaffer (1982) tested for 
displacement forecasts. Shaffer (1982) and Shaffer and 
Elsberry (1982) showed that for tīnis inispendent sample of 


meeeica: Cyclones, the forecast error from a technique based 
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TABLETI 
EE o e e a a a 500 mb. (shatters 1962) 
Mode. Eigenvalue Explained Variance 
1 39.86 33 
2 21.60 51 
3 9229 59 
y 7.43 65 
5 6 .08 70 
6 529 75 
7 4.00 78 
8 3.13 81 
9 2-49 8 3 
10 2 4 85 
15 1:203 9 1 
20 „ale 94 
40 .11 98 
60 204 99 
129 «90 100 
on an EOF representation of Sie ASADO... ) fOesing is 


decreased by an average of 175 “hen compared to the official 


JTWC (Guen) forecast error. 
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15 20 
Begurce 3. Eigenvalues of the 500 mb D-value Fields 
(CASES Mana” cae done carlo s35neraved  . 
Ezgenvalues (“c21angli2s) Based >n 100 Simuiatiors 
oí Random Data. 
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III. DEPINITION QE RECURVATURE 





The definition cf a recurvaturs case can be thought of 
as a Set of restrictions based on the observed track of the 
Bropical cyclone. The definition of racurvature that is 
used in this research is that described by Leftwich (1978, 
1982): 


"Recurvature is a net displacement northward of 3159 
during che forecasw period al the ang of 
DemensastwatdemMotmanm during che 12 kours prior to the end 
of tha fcrecast period". 


M mns are net considered if thair past 12 hour’motian is tc 
the east. A graphical dəscriptrion cof  recurvature and 
non-recurvature cases is given in Fig. 4. 

Some cf the cases present rather lifficult forecast 
Berne, Figure 5 denicts Typhoon Trix (initial pcsition 
E )o41T 27 August 1971), which for tha initial 36-h period 
would bs considered +o have recurved, whereas for the 54- 
and 72-h period it would be considered as a non-recurvature 
case. PaeacmomGeprces Typhoon Polly (iritial position at 
12G4T 29 August, 1974), which for “che 236=-h period is 
considered to have  zecurved, but then moves westerly to a 


DS position which 3s considered as having not recurved, 
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LEGENO 
4 * NRCV 
° = RCRV] 
° = RCRV2 
+ - NOTC 


25° N 


=z 
e 
110° E 120° E 130° £ 
Figure 4. Examples of Ideaiizel 
recurvature (3) and a 
and finally moves north-northw 


recurvature case at 72h. 


fen, 02 October 1975) shown in 


racurvature case at ail three 


nuce (initial position at 12GMT 


8, is a non-recurvature cease for 
Anita (shown in Fig. 9) woul: 
12GMT 22 July to 00GM1 23 July 


EE ward notion of the tropic 


considered nor The 12GMT 23 July 


A summary of the 


storm sample is qiven in Tabie 


24 





recurvature 


LM e 


140” E 150° £ 
Recurvatuze (1,2), Non 
Not Considered Case (4) 
est to be considered as a 
INS cOla “(on2 claw position of 


Hag. 7 is an ¢xampie of a 
periods, whereas Typhoon 
16 September 1975) in Fig. 
all three periods. Typhoon 

not bs considered for the 
iaa posicions dus tc the 
aL CyCo, Duc" woulgdg be 


to QO0GHT 24 July period. 


Statistics for the 


IL. The 
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Figure 7. Typhoon Cora Case Illustrating a Common 
Recurvature Case 


LEGEND 
SiS 
+ = ALICE BEST TRACK 


laa lap aaa aja 3 


4 
0 
à 
r 
4 
b 
0 
1 
0 . 
4 
t 
1 
0 
4 





110° € 


rigure 8. Typhoon Alice Case Illustrating a Non-recurvaturs 
case 
statistics were obtained in the same nanner as defined 


earlier, except that the motion during the past 6 h was used 
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115°E — 120'E 125* E 130° E 135° E 140° E 145” E 
gure 9. Typhoon Anita Cass MR cad the Exclusion of 
the 12GMT 22 July to 0QUGMT 23 July Periods Due to 
the Eastward Motion 


for the second consideration for recurvature vice the past 


teh motion. All past meridional and zonal motions were 


tt 


@eeained from the warning positions 9 sach individual 


te actual operational 


Ul 
pa 
tj 
a 
f» 


storm. This was done to 
Bemazttions (Shaffer, 1982). 

À graph cf the recurvature statistics showr in Taole IT 
is presented as Fig. 10. There is a slow increase of 


probability of recurvature with tine. This increase may be 


(0 


E 2sesd by the follcwing rsasoning. A tropical cyclone in 
ts developmental stages is mor2 often in the easterly rađe 


wind regime and moving west of northwest (non-recurver as 





TABLE IT 


E A EEES EEC 
Recurvature (prob) for the 504 Desi eal Cyclone Sample 
5... ....- tau rcv  nrcv  nohist notcon prob = 

06 131 234 0 139 HYS) 

12 125 240 0 139 MEI, 

18 130 239 0 139 338 

24 134 231 0 139 S sk 

30 135 230 0 139 23709 

36 158 196 22 128 .446 

42 119 169 115 101 2413 

48 125 163 115 101 .434 

54 133 154 1 16 101 .463 

60 133 139 142 90 .489 

66 96 111 290 67 464 

E 103 104 230 67 .498 

75 107 100 230 67 E 

34 99 87 254 64 215132 
Pez ned) . MS theftrop cal cyclone 1nt=ensifies and moves to 
higher latitudes, ae begias to interact wich the 


mei latirude Westeriies and therefore is more likely to 
recurve. Cheng-Lan and Sadiar (1982) describe tihe mes 
frequent track of a recurving Tropical cyclone as "a Tegular 
Clockwise parabola". That is, as the *ropical cycicne moves 


WeStwWard it also moves northward, with 2 greeter likelihood 


EE cung north of a 3459 lina from the initial position. 
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PROBABILITY OF RECURVATURE 
0.00.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 


a a.) La E T OT AT nÀ Á 
6.0 12.0 18.9 24.0 30.0 36.0 42.3 48.0 54.0 60.0 66.0 72.0 78.0 84.0 
TIME IN HOURS 


Figure 10. Probability of Recurvature (504 Case Sampie 
Ciimatolo gy) 
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The statistics contained in Table II and Fig. 
designared as “he sample climatological pronabilities oí 
recurveture. I= should be mentioned that the independent 
Sample statistics for  recurvature for all time periods are 
Within 7% of the sample average probability, which indicates 
consistency between the two sanples. Ihese values wiil be 


- 


used later as a standard mor =. valieron 5t varsous'forecsss 


Schemes in predice ing  recurvatuse at specified time 
Enteryals. That is, a scheme will be said to have skill if 
lE predicts recur vature more accurately shan tpe 


meeemastOlogicel probability. 
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IV. REGRESSION ANALYSIS 


Neumann and Lawrence (1975) demonstrated that a 
Sta istical-dynamical modal using prognostic and 
persistence-related variables in a nonlinear regression 
scheme showed much promise in forecasting Atlanti ocean 
Meopical cyclone motion. Shaffer (1932) also used 2 
statistical-dynamical scheme that used 10 BOF coefficients 
to represent the geopctential field (analysis time only) and 
a linear regression scheme to predic eropreal cyclone 
Met cn. The SOF analysis Hot only sänpli-ies the 
representation of the synoptic forcirg, it also reduces 
Enun»t2cant the cost of producing an accurate and viabie 
forecast. 

en npor ance of the nonlinear terms in regression 
scnemes has been demonstrated by Neumann and Lawrence (1975) 
and Paegie and Haslam (1982). Although nonlinear +erms were 
included in those studies, eeanomleseeons:deratsons and the 
ENUSIJCc2TY Of application of a linear model led +o its 


eaopc20n in this study. 





DNE cs ne (effect Joumeate Luding EJF coefficients as 
pređictors of tropical Cyclone motion, three E er oals 
for each of the prescribed times were completed. The first 
test (refered to as PERS) was accomplished using only the 
persistence-related parameters and the initial latitude ane 
longitude. Ac deseruapeMon of ali potential predictors is 
contained in Table Ill. (All laexymma szene longitude 
displacements (displ) are in nautical miles.) The second 


experiment (referred to as EOFS) utilized 10 zOF 


(D 


eef ficients and the initial latitude and longitude. Th 


esi Sctence and 


third experiment used a combination of the p 


(D 
D 


EOF predicters and ís noted as P&E. Tt should be mentioned 
that the cyclone maximum sust2ined wind svesds at 12, 24 and 
Emmet DELOL tc the initial observation wsrs also incorporatad 


as possible predictors. This was dons x» accordance wit! 


rn 


Riehl's hypothesis (1971) 2 maximum intensity being 


Mmasocre-cd With zecurvazure for Atlantic tropical cyclones. 


D 


However, none of these additionai predictors were chosen 
during the stepwise regression. 

The 504 case tropical cycione sampl2 was analyzed for 
recurvature as descrired in Chapter III. A value of 1.0 was 


assigned to the storms that recurved and a value of 0.0 was 


(n 
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ene ou ee 


Potential Predictcrs Used to Da AS Regression 
Equations for Each of the Three fime Intérvals 


Potential Predictor Name Description 
Variable Number (Coeficient-Coef) 
1 COF1 Eigenvector 1 
. Coefficient 
2 (e sigenvector 2 
„. Coeerıeyent 
3 COF3 Eigenvector 3 
~ COS = tue went 
- COFU B GSV C Or 4 
SEDer Tr 2eIent 
D COTS Szgenvseeor 2 
OS EU nt 
6 COTO Elgenvecior 6 
. SOSELa cient 
7 CONDE Eigenvector 7 
~- Coie mime Lent 
8 era Eigenvector 8 
eon LG Lene 
B Goro rigenvector 9 
, COS ee ot 
10 (Om Ez:cenvsee or. 10 
$ ICAA eee 
11 PLATO AE O oe 
Latitude 
12 PLATI Past 12 hour 
meridional displ 
13 PLATZ ‚Past 24 hour 
nezi.dr Snalzeavepn! 
14 PLAT3 Pease sso hour 
meridional displ 
15 PLATY Meridionai displ 
Or pa loc 2=-2U irs 
16 PLATS ecca onana 
foc pT or 12-30 0TrS 
17 PLATO Merlo al a Sop. 
COLE OC 2U=36 MTS 
18 PLONO nus a LO Cone 
lcnqitude 
JE PLONI pora Eu 
AN pL 
20 PLONZ2 Past 24 hcur 
zonal displ 
21 PLON3 PUSO. MOE 
zotal dis 
227 PLONG Ona leeds OL 
fO prior 24 hrs 
23 PLONS Zonal displ 
BOPREDLITOP 12236 hrs 
24 PLON Zonal displ 
fcr prior 24=36 hrs 


assSianed to the cyclones that had not recurved at the time 


Of interest.  À value of 999. was assigned to cyclones that 


MERO history at the time period of interest. 
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The BMDP stepwise linear regression program screened the 
dependent data set for storms with past 12 h longitude 
displacement (PLON1) that was to the east and also tropical 
storms that had no history. This screening accounts for the 
different number of cases considered for the regression 
routine and the climatology cases at the specific time 
periods of interest (see Table II). The equations at 36, 54 
and 72 h (see Tapies IV, V, VI below) are derived from a 
stepwise linear regression using an F - ratio of 4.0. The 
F-ratio is a test of the significance 55 the independent 
variable coefficients in the regression squa*ions. The 
number of cases of the dependent sample that are selected 
EE cach regression period is noted in the ccorrespondin 


table. The number of cases decreases with time due tc the 


df 


Mas trictions or inclusion discussed in Chapter II. Of not 


har most of the cases at 72 h are 3 subset of the 54 kh 
Sema Set, which is in turn a subset of the 36 h data set. 


iD 


MS E Squared Statistic ,RSQ, (percent of variance in +h 


pa 
o 
D 


predictarG explained by the ragression =quation) wil 


Meea for discussion of the results. 
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The linear regression equations 


potential predictors of recurvature or 


SE the rorm: 


P(0:1) = a0 + a1X1 + a2X2 + a3X3 + 


Tables IV, V, and VI list the coefficients 


for each of the three tests generated from 
data sets. The number of cases used 
coefficients are 259, 206 and 138 respectively. 
TABLz IV 
Regression Equations for 36 h Recurvaturs = 
Three Methods Tested (Dependent Sample), an 
and the Standard Deviation from the Mean (STD) 
PERS DUNS PS E MEAN 
intercept .1852 =sco/] ASI 
PLATO .0167 19.4 
PLONO .0054 13322 
PLATZ 20017 . 0014 96.4 
PLON1 0031 0031 -92.6 
COF 1 -.0326 2:02 15 . 81 
COF 3 200357 -.24 
COF 5 -. 0 302 2520) 
COF 7 20341 02 99 0 
COF B .0493 =.) 
COF 10 -, 0545 -.0397 -. 10 
Bach regressicn coefficient is unique =xcep* for 
Of PLON1 fcr the PERS and P&E test at 36 hours. 
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produced 


and 


to generate 


from tne 


nor-=recurvature are 


* a7X7. 


variables 


she dependent 


the 


or Each o£ the 
^ Mean Values 


Dre ark) 
1223 (0%) 
90.2 (NH) 
55.8 (NM) 
Sez 
2.89 
DIEM 
1.91 
w77 
1.48 


the value 


pto hou" 





TABLE V 
Regression Equations for 54 h Recurvature as Described in 
Table IV 


4-— ap a D op UH AG A umm m —D EP aD «ND wb GEB enm EDO ED amm y AO A A AUD EB e UE END A SS FPSB SPP See eet A aA umb eSB Be A Ss Ss =F m =» 


PERS GOES Eo MEAN STD 
INTERCEPT .4446 ©4957 24873 
PLAT2 . 0019 .0014 J976 86.6 
PLON 1 «0022 . 00 16 zs 26.6 
Gor 1 2200312 220252 lcs 4.8 
COF 3 03877 20737 290913 2.74 
Cer 8 .0618 .0465 E 1.32 
COF 10 -.0514 US 1.50 

TABLE VI 


Regression Equations for 72 h Recurvaturs as Described in 
Tapl2 IV 


< da A 2 A APA A A OPS AED A A UND ED A SPs SP ewe BSP TP A A A A Ay AA A SP Be see ses BP SEE y A A A <i> SU: Bs SEM ee PE Se = 


PERS LOPS PS E MEAN SL) 
intercept «5487 5398 25921 
PLATU .0028 200.26 44. 1 48. 1 
PLON 1 0027 . 0024 -89.1 222 
COF 1 -.0324 2902090 2209 055 
EOF 3 .04 10 =. 75 LEY 


be noticed that the value of PLONI must necessarily be less 


un 


than Zero, because a PLON1 greater than zero represen: 
easterly movsment, which is considered as having already 
recurved. 

The totali number of combinations of predictors and 


regression coefficients wiil not be discussed, but several 


combinations will be mentioned. EE (pase 24 h atitude 
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displacement) is almost always greater than zero (the 
exception being trochoidal oscillating storms), and makes a 
considerable contribution towards Tecurvature (pradictand 
close to 1.0) especially for large displacements. PLATO and 


PLONO {initial latitude and longitude) are always greater 


Ban zero, witha larger values being favorable for 
recurvature prediction. Piece, COCMETel ent Mr the» first 
eigenvector (COF 1) has the opposite effec. Physically, 


the stronger the zonal pressure gradient (large COF 1), the 
lower the probability of recurvature. The opposite is also 
mue; the prcbability of recurvaturs is higher for 
extremely weak zonal pressure gradients as illustrated in 
pug. 11. 

Analysis of the sample mean values of the coefficients 
SEO 1 (assuming a normai distribution about tne mean) 
catas that the number of  recurving tropical cyclones 
would be less than che number of non-recurving storms. POF 
example, the number cf recurving storms at 36 h (T2blé IT) 
is 158 and the number of non-reacurving storms is 196. The 


Seetticient of £0F 5 is also quite interesting. The more 


ct 


the EOF 5 pattern describes the synoptic forcing (the larger 


the value of COF 5), Pie mocro likely the storm willi nox 
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pure 11. Storm Tzacks Relative to the Initial Position 
Zora) oeeets wien an Orthogonal Coel icien: 1 
tha lS (a) Between -9 and -39; or 52) Between 
Pomona Pech CFOSs Indicates a 12% 
Displacement (Shafier, 1982) 


recurve. The converse is also true in that the stronger the 
inverse or tne negative of tne pattern (a strong shortwave 
between two ridges), cierro oe kely the cyclone is =0 
recurve. 


Ms XS»NUO" usans a complete discussion of ali th 


(D 


Peesciois 54 linear combinateons of predictors, coefficient 


N 


and their inverses. It does show how a simple combination 





of EOF patterns can be applied with straight-forward 


meteorological reasoning. 


A. 36 H REGRESSION CASE ANALYSIS 

As expected, the predictors selected for the PERS 
regression analysis included the past 12-h longitude 
displacement, the past 24-h latitude displacement and «ha 
initial latitude. Simple extrapolaticn techniques do well 
for these relatively short-time forecasts. The longitude 
displacement over the past 12 hours raduces the variance in 
the predictand (increases RSQ) slightly more than the 
latituda displacement. The selection of initial latitude as 
HENENEUTeCdietor is consistent with the vVarliatiom of the 
probability of recurvature with time (and thererore latitudes 
for an idealized storm track) as discussed in Chapter Il. 
The total explained variance of the predictand for the 
persistence-related predictors «as 29.2%. 


The individual EOF predictors reduce less of the 


(1 


(D 


variance than the persistence-type predictors (Table VII). 


cr 
me 
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Ce that eigenvector 1 (Which is an indicator of th 


strength of the zonal pressure gradient and large global 


rth 


wave number} causes the largest reduction in the variance o 


the predictand. The remaining five eigenvectors that wer: 


SE 





selected are indicative of smaller scale forcing f2atures, 
and reduce the variance of the predictand much less that ine 
first eigenvector. Shaffer (1982) indicated the possibility 
of weak northwest motion (or southeast motion for the 
inverse mode) would ke associated with eigenvector 7 anda 
strong neriodional forcing being associatsd with eigenvector 
Bre Of special interest is 2igenvector 5. Although the 
reduction in variance associated with this eigenvector is 
small, its inverse pattern m Ssinilar tö a strong short 
wave in the westerlies which could advect the storm 
northward. The last predictor selected, eher initial 
longitude, again is consistent with the variation of the 


recurvature probability with time. Ihe “cotal” or the 


ral 


e lained variance for the cO0rS test was 26.5% , which is 
2.7% less than using persistence only predictors. 


for tha Žinal test Shows à 


Ui 


The selection of vredictor 


combination ci both PERS and EOFS predictors. The largest 


O 
I 


Beauction of the Variance occurs with past noti 
(persistence related) parameters. The strength of the zonal 
pressure gradient appears as ponerte most importan: 


features (EOFS 10 andá 7) 


(1 


puedictor with smaller scal 


contributing smaii increases to the explained variance. The 





total of the explained variance of the predictand is 34.0% 
which is greater than either the PERS or EOFS methods, but 
Perl only a small fraction of che total variation. 


Part of the reason tha: more of the predidictand 


variation could not be explained may be related to the 
wo ar y definition of recurvature, 2.9. a movement =o 
320° is recurvature (value 1.9) whereas 310° is not (value 


of 0.0). When a yes/no decision is applied across a line as 
in this situation, one cannot expect an explained variance 
near 100%. The me*hcd could stili provide ussfui guidance 


puce orf Salz. en rozecaste is 
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J 


ror the forecaster if the perce 


higher than other metnods availa bie. 
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The entry of PLATO, and then the subsequent deletion 


esz3on oF 


qt 
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| se 


nd 


(25 
fdo 


ee OF information is added, is an 
overlapping information. Notice (Table VII) that PLATO dees 


rot enter in scheme (c) at all. 
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Maps YI 


Summary of Stepwise Linear Regression Results for 36 h (259 
cases) R Squared (R50), Increase -n R Squared (RSQ I ) 
and P to Enter or Remove 


P.I > a wb ab wb qb qb q qb qb qD wb es A Y A A A OE A» aes BSP ese Sas Se Fa SEM E E ER JS see E O SPS SP y y eee A A ee ADA E o 


a) PERSistence only 


Step Parameter ESQ Hoo Iner F enter F remove 
1 PLON 1 151 «5 1 45.65 

2 PLAT2 «120513 :1 12 332.06 

3 PLATO 27292 4/0129 10.42 

b) EOFS only 

1 Com «21 2121 3200 

2 PLATO . lo .040 12.04 

3 COFEB ao Y 26 8.20 

4 COr7 . 206 «Uo 6.04 

> GOF3 s227 219241 6,85 

6 COFITO . 243 . 0 16 sy, 

E PLAT e 234 -. 009 2.94 
8 COF5 . 288 .9 14 4.54 

B PLONO . 205 917 Se 5 

C) Persistence & EOFS 

1 PLON1 sil „is 45.65 

2 PLAT? + 203 «d 12 39.06 

3 COF] «318 2090 18.54 

4 CCF10 az ~914 De YU 

5 COF7 . 340 sah 9a 18 


Be 54 a REGRZSSION CASZ ANALYSIS 
The selecłŁżion of EERS predictors at 54 h {Table VII) are 


Ee same as e- 36 h (Table VI) except Zar the deletion of 
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mr initial latitude. Had tha F-test been 
Ino. the initial  1iatitude would have been  seiectsad 
subsequent to PLONÍ!, as was tha case for 36:h. At 54 à the 
latitudinal displacement for the vast 24 h decrsassd the 
variance of the predictand more than did the vast 12 h 


longitude predictor. These two predictors comprise 16.4% cf 


4 1 





the variance of the predictand, which is almost half the 
value obtained at 36 h. This decreas is not surprising, 
because simple extrapolation tachniques have less validity 
Bes 4 than at 36 h. 

The number of predictors selected in the EOFS test 


decreased to four, ccmpared to seven at 36 h. These four 


meediceors at 54 h are all included within the first fiwe 


E-dic-ors a- 36 h. The north/south pressure gradient (EOF 
NN again contributed to tha largest reduction in the 


Variance of the predictand (see Fiaure 2). EOF 8 a 


KD 
"O 
(D 
W 
13 
{N 
ct 
O 


heme Similiar to ZOP 5 in that it is assocsated with a 


mid-latitude short wave combined with an area of high 


r 


pssucze dirscily south of the tropical storm. EOF 8 is 
Seereiated airectly wach recurvature as are EOFS 1 and 3. 


Mere 10 LS andirectiy correlated with recurvature due to tne 
PA pressure area to the northwest of the storm center. 
The stronger this pattern, ths less likely it is that the 
storm will recurve. moe tota! dort the sexpleined vertance of 
the predictand is 19.4% which is smaller than the value for 
DO using EOFS only, but it is larger than 54 bh regression 


using persistence parameters only. This result suggests 


that synoptic parameters begin to explain more of the motion 
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at later times than does the extrapolation associated with 
persistence variables. 

The combination of the two types of predictors expiain 
24.6% of the variance of the predictand. The only EOF 
predictor that does not appear in P&E is eigenvector 10. 
The inclusion of the eigenvectors as predictors increased 
the RSQ amount by 8.2%. The decrease from 34.0% expiained 
variance at 36 h to 24.6% at 54 h was expected due to the 
likelihood of more tropical cyclone interactions with tne 


synoptic environment. 


e 


PRBLE VIII 
Summary of Stepwise Linear R2gression Results for 54 h 
(206 cases) as Described in Tabie VII. 
a) PERSistence only 
Step Parameter ESQ RSO INOT F enter F remove 
1 PLAT2 - 102 « 102 2329205 
2 DLON! . 164 . 062 15.06 
b) EOFS only 
1 COF1 e C74 <. 0 74 16.20 
2 COF ol Vee 20:9 3 12.42 
3 COF3 - 170 043 10.41 
4 COFI0 . 194 .024 6:93 
C) PersiStience 6 ZOPS 
1 PLAT2 102 492 23.226 
2 PLON1 . 164 1, 62 13-056 
3 COFI «203 USO ee 
4 COFEB > 230 “021 Cea? 
3 IT . 246 NS 4.15 
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GF 72 KH REGRESSION CASE ANALYSIS 
UNS Selección ot predictors for the 72 h prediction 
includes a new predictor, PLAT4Y (latitude displacement fcr 


mre period 12 to 24 h prior to the initial time). Lass 


(D 


predictor is inherently positive and therefore adds to th 
probability of recurvature. The past 12 h Ilengituds 
displacement is the rémaining predictor chosen. Tue oval 
of the expiained variance using these two predictors is 
13,0%. Although this is less than half of the value of the 
BH regression, it is only slightly less than the 16.5% 


obtained for the 54 h regression. 


TABLE IX 
Summary of Stepwise Linear ee ression ne Tommi? -h 
(138 cases) as Described in Tabie VII. 
a) PERSisten only 
Step Parameter ESQ RSO nes F enter F remove 
1 PLON 1 40/28 2:970 10.84 
2 2LATA «2159 «0 76 12291 
b) FORS ONLY 
1 COF1 . 076 . 9 76 1505527 
2 COF3 1238 2310597 203 
C) Persistence & EOFS 
1 COF1 "00% O Des), 
2 PLON1 zer 056 8.63 
3 PLATY . 192 «0:60 9.94 


Tne EOFS-enly regression produced results that explain 


moon Or the variance of the predictand than dia the 
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persistance-type predictors. Coefficients 1 and 3 explain 
12.8% of the total variance of the predictand, Week Ss 
almost identical to the value for the first steps of the 54 
h regression. dowever, two different coefricients entered 
the 54 h regression equation. 

The combination cf two groups of predictors is again 
superior to either ipday i dual group. Wee 12 h the 
north/south pressure gradient is more important than the 
longitudinal displacements 84 h earlier  (PLON1) or the 
Meettude dispiacement 84 +59 396 h prior to the initial 
Position (PLAT4). Use of extrapolation-type methods can not 
be expected to contribute as much as the synoptic parameters 


at such long time intervals. AS the tropicai cyclore moves 


4) 


comparatively less 


Ui 


anto the mid-latitudes, persistence i 
Userul as a predictor than in the tropics. 
The amount of expiained variance of the tecurvature 


event decreased with time for 22ch of the three test groups. 


the predictana) is the ccmbination of parsistence-related 


(1 
an 
iD 


alles ana the EOF coefficients which describe 
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ct 
ps 
bet 


synoptic patterns. Although the PSE technique consisten: 
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ce 

3 

O 

Qs 

Ul 
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Explained more of the variance than the other two m 





further testing is necessary +o determine if the EOF 
predictions are superior to other techniques such as 


climatology or mere chance. 
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V. TESTING THE REGRESSION ANALYSIS 


Choosing a forecast verification scheme is very 
important for testing the ability of a regression model or a 


Peognostic method tc forecast a specific event. The 


tty 


geophysics community is just one of the many scientific 


> 


communities that desire an accurate and acceptabie method 
for describing the skill of a given model. Muller (19443) 
provided a literature survey of 55 short-range weather 
verification methods that were developed during 1884 to 1943 
in six different countries. Muller divided thea schemes ir-o 
two areas: comparisons of the forecasts and observations as 
a ratio (percentage ccrrect) and verification methods which 
w compared to chance or climatclogy, the so-called ‘biind' 


forecasts. 


4) 
O 
ti 


In the early verification studies, “hs schemes 


rh 
O 
17 


comparing forecasts with observations had +o be simple 


mime case Of applicatilcn. Verification methods today compare 


BeEDredict2zon skill against some clinatological value or to 


t 


pure chance as a point of reference. 
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Three skill scores will be iiscussed: she Heidke Skill 


th 


Ecore (ASS), the Brier P=-Score (BPS), and the percent o 


correct forecasts (PCF). The Heidke skill score can be 


written as: 


Hoo = Vee Cc 





where VF is the number of correct or verifying forecasts of 


an event (either correctly hforacast recurvers (CFR), i.s., 
predictand is greater than or equal to an assigned 


anma ology probability (Table IT); or a correctly forecast 


non-recurver (CFNR), i.e.,„ the predictand is less than the 


un 


Em mMatology probability), T is the total number of forecast 
made ard C is equal tc the number of correct forecasts that 
would be expected from chance. In this test C is always 


one-half of the total number of forecasts. Ties woule 


) 


jd 


t 


have a maximum value of 1.0 for perfect forecasts 4 
value or zero when the number of correct forecasts is equal 
[ohne number of forecasts expected to be correct by using 
Chance. Negative values occur if the nethod is inferior to 


the comparative methods. 
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The Brier or P-Score has been widely accepted as a 
method to test probability forecasts. The binary form of 
the EPS can be written as: 


N 
BPS = 2 * >) Fest (4) - ACT(i)) **2 
NE oe, 


where FCST is the regressicn equation predictand and ACT is 
either 1.0 if the storm actually recurved or zero if 
recurvature did not occur. A perfect BPS is a value of zero 
and the worst possible forecast wculd be indicated as 2.0. 
It is readily apparent that the BPS is actually the nean 
Square error cf the set of probabilistic forecasts. 

The percent of ccrrectly forecast events (PCF) can be 


Weatten as: 


PCPS CER + CENR 
T 
where CFR, CINR and T are the same variables as described 
earlier. 


As discussed in Chapter IV, the preälctanä values should 
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lie cn the continucus interval betwe 
unfiltered predictand statistics for both the dependent and 


independent samples for all periods appear in Table IXI. 
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Values larger than 1.0 (indicative of a strong recurvatur 
Signal) and smaller than 0.0 (Erdredeing strong 
non-recurvature signal) exist as maximum and minimum values. 
These values outside of the defined range are due to the 
assumed linear relationship between the predictors. Table 
XII provides similar statistics for the  predictands after 
values that were more than 1.0 and less than 2.0 were ser 
equal to 1.0 and negative values were sat equal to 0.0. 
Comparison of the dependent and independent samples of both 
tables indicate that the two s2ts are quite similar. The 
different number of cases included in several tests is due 
to the predictors PLAT2 and PLON1 being =xcessiveliy large, 
(purposely assigned flagged values) which causes several 
predictanā values for PERS and P&E to be in excess of 2.0, 
me upper filter limit. 


EBENEN C Yy»099nd Xwindicat2o tha: thə variance or ths 
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dependent predictand is great 
variance of +he independent  predictands ana 15, of T8 
eompar Sons. This is evidently due to the small size cf the 
independent sanple. 

It also should te noticed ‘that +*+he dependent sample 


Variances are largest at 36 h and decraass time. Tn all 


(^ 
tJ 


50 





dependent data set comparisons, the variance of the P&E test 
is larger than the other two tests. This is most likely due 
to the 'sigral' of either rscurvature or non-recurvature 
being strongest for this set which includes both persistence 
and synoptic features. The smaller values of the variance 
me fable XII, as compared to Table XI, is due to the 
filtering of the predictands, which is described above. 
TABLE XI 
riltefing. The Number of cases, Hein, mikinuh and Mininun 


Values (mex,min) and the Variance (var) of the Predicted 
Values Mor fhe Dependent and Independent Data Sets are Shown 


Dependent Set Independent Se* 
PERS EOFS PES. DERS EOFS PEZ 
cases 349 433 350 48 ug 48 
36 h mean (i98. (433) (233! Gel O36 Cis) 
36 h max 121907 14. 550 810 1.020 7732299 ,,12059 
Bo A min -. 256 -. 243 -.259 =a 299 ed] 
36 h var e. 120 . 083 «153 .096 42063 084 
cases 271 345 271 4 3 43 43 
54 kh mean (ide) (32) rene Gb a3 (iz 
54 h max ID IEA 12790 SE Ec E 
seh min = 2021 + GI —, 112 = 92 =.,392 
5S4 h var -065 . (49 .990 4.00 7 ou SONS 
cases 176 235 176 38 38 38 
SES wan 885 $27 Qs ss Ga 6» 
72 h max 1.360 1.06060 1.349 .978 .821 TU ^U 
mh un -.037 . 099  -.038 -.048 see = SZ 
M2 h var S053 202 20:65 053 . 028 2062 
The summaries cf the HSS which compares the threes 
Bemeds oz *esting using the climatolsgy probabilities of 


recurvature for both the dependant and iniependent data sets 


are presented in Table XIII. LE snould be noticed that the 





TABLE KILL 


Pee dietand Statistics for 36,.54 and 72 h After Filtering as 
Described'in Table XI. 


ED A E p up AA E EY A ap ap ap HEP ED ap SPS See ee See D Amp p SUED XD eee eS eee GM SSP SS SPSS SSS SS Se SSeS SS SS SS A TO A e SEE oe 


Dependent Set Independent Set 
PERS pe PEE PERS BOFS P&E 
cases (349) (433) (350) (48) (49) (48) 
36 h mean . 486 . 427 . 4 38 . 466 . 430 ung 
36 h max ie 000 1.000 12000 1.000 1.000 1.009 
36 h min . 000 . 000 .000 . 000 025 -000 
36 h var .088 2069 «099 -082 2 054 20075 
271 345 274 43 43 43 
EE mean (i29 (333) 2) Cras Crab Gar 
354 h max 1.000 1.009 1.000 4867 4983 -3889 
54 h min . 000 . 000 .000 . 000 $009 2900 
54 h var $051 . 046 .068 . 042 +. 939 ES 
cases 176 233 175 38 38 ES 
72 h mean 33? (37 (43) 6) al Gi 
72 h max "COQUE. 000 1.000 -978 2021 1.000 
72 h min .000 un .000 . 000 a12] +900 
72 h var .049 027 2057 -052 . 028 2:057 
number of cases in each test does not change from the 
earlier values. 
TABLE LILI 
Heidke Skiil Score Usina thz Sampie Climatology as a 
Sedad CE Compa? isor for Defining Recurvature Versus 
usu-rscurvature for the Filtered Prediciands (as in Table 11) 
ARESO OCET O TSwWenzecrn, 
Deperdent Set indeverden=- Sex 

PERS BOFS PGE PERS EOFS P&E 
36 h >53 ol, 2926 .464 -296 . 469 
54 h e505 2 422 1965 .442 2536 e440 
ES h .431 . 254 e454 2431 213163 5.10 


NUurosnd-ion . recte-p*ng to the skill of the HSS it 


I 


Table XIII must be carefully interpreted. When using th 


NS. any predictand that is greater “han or equal to a 
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specific value (in this case a climatological percentage of 
recurvature for that forecast interval) is given a forecast 
recurvature value of 1.0. Any predictani that is less than 
the climatological value is considered O be a 
non-recurvature case, and is assigned a value of 0.0. These 
'newly assigned? predictands are then compared +o the actual 
storm zrecurvature values. It the forecast and the actual 
value match, a verifying forecast is tailied (CFR or CFNR is 
incremented by one). 


A graphical comparison of HSS versus a varying 


Fh 
O 
ry 


climatological value cf probapility of recurvature (used 


(1 
Ly: 
(D 


comparison) is shown in Fig. 12. The maximum values in 


(i 


ESS are close to, although not necessarily associated with, 


meme climatological probebility values indicated for the 


Q 


three different time periods. 
Another point needs to be mentioned concerning the HSS. 


twee climatological prcbability values ror zecurvature (Table 


II) at the specified times should be conpared to the means 
and variances of the predictands shown in Table XII. Al. 


three climatciogical values of recurvatura fall within the 


ABS? standard deviation around the mear in all three tests. 


1 


BEEstnis An Mind, ore would oredict tha: mechods with the 
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LEGEND 
Q= 36 H P&E 
© - 54 H P&E 
4-72 H P&E 


HEIDKE SCORES DEPENDENT SET 


-0.20.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.6 0.9 1.0 
A ee 





0.050.125 0.200 0.275 0.350 0.425 0.500 0.575 0.650 0.725 0.800 0.875 0.550 
PROBRBILITY OF RECURVATURE 


Figure 12.  Heidke Skili Scores with the Value of thea 
Drew -dnmgErronerhs Poe Equacions that is 
Pages” tc be Recurvature Varied from 0.05 to 


smallest variances 
greater probability 


(low signai regime] 


would have a 
orbe ng n2ar ne 


than 15 


higher 


HSS due” zo. *he 


Climatological value 


*he tails (high signal regime) 


Ehe predictand distribution (to be discussed later). 


En Table XIII, 


"he combined P&E 


highest scores relative +o climatology, 


period with 


best. An 


-he mean appears 


tne dependent set, 


Mets On voz 


where 
she clustering 
the 


54 h 


a 


Ea nd 


54 


method produced +he 


the 36 h 


in 
Ane 
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method was 


about 
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predictands based only cn EOFS have the highest scores 
relative to climatology at this time, while at both 36h and 
72 h the EOFS were infericr to the other two methods. P&E 
generally showed the highest skill for the independent sex, 
with the exception being at 54 h aS discussed above. 


The distribution and cumulative frequency of the 


(D 


filtered predictands for the three tests ars given in Table 
IU. XV, and XVI. An unbiased distribution would have 10% 
of the observations in each category. This ideali lima: is 


nearly attained at the 36 h time period for the dependent 


data set. However, biases become more evident for later 


taj 


time periods. The FSZ test genetaily is less biased than 
the PERS test. Both P&E and PERS are less biased than the 


Sites, which again is an indication of the clustering of 


the EOF predictands néar the mean value. 





TASES KIV 
eee soatcrand DiSt= ibu<-ion and au, cn tage 
Meeri ibution for 36 n. Cumulative Percentages 2 Ba rens. 
Dependent Sect Independent et 
B6 h PZRS EOFS P&E PIRS BOFS P&E 
> «| 36 (10 38 (9) 41(12 BOT 3 (6) 52/010 
E .2 37 (Z1 BUE 31(21 OAT 5 (10 Sl 
Ec .3 34 (31 70 (37 4233 SN 6 3) 4 (25 
Ea .4 35 (41 69 {53 40 (44) 8 (44) 6 (41 od 35 
2 .5 30455 47 (66 42 (56 à (52 12 (65 9 (54 
E .6 32 (64 43 (76 31 23) 8 H 8 (82 8 71) 
Eo .7 33 (7ü 33 {383 23 (71 4(77 5{92 5(31 
> .8 LOS 26 (29 22 (78) 3 (8> 0(92 3(88 
.8- .9 14 [30 16 (93 23 (84) 6 38) 1(94 1(90 
A 1.0 44 (100) 31(109) 55(100) 2(1090) 3 (100) 3 H5 ) 
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TABLE XV 


Predictand Distribution and Cumulative Percent Sale 
Distribution for 54 h. Cumulative Percentage is in Parens. 





Dependent Set Independent Set 
54 h PERS EOFS P&E PERS EOFS DEE 
mo .1 17073 4 1 6 4{ 9 2 
n- D2 20(12 4 13 22 14 2(14 1 3 3 18 
= 3 30 (23 59 (30 37 (28 3424 6 (21 gl 19 
moe « 4 45 (39 61( 47 42 (43 539 8 (40 T’ 
Ni .»5 43 (55 63 (66 42 (59 14 (65 10(63 6 (49 
Bo .6 39 (69 51(80 act T2) 6 (73 6 (77 D 
Eum .7 39529 2 26 (€8 21 (82 6 33) 6 (91 Sis! 
1 .8 20,09% 19 (93 14 2 1(95 3 28 2195 
Eu .9 13(97) 11( 97 9119 51199) au 22110) 
n9-1.0 9(100) 12(100) 28(190) 0(109 1(100) 2 (100) 
TABLE XVI 
[en Distr turen and Cumulative Petcentage 
Peeseribution £ 2 n: Cumulative Percentage is in Parens. 
Dependent Set Independent Sat 
72h PERS SOES PSE PERS LOGS P&E 
Nu- . |! 6( 3 1{.4 11( 6 3 085 0( 0 2C 09 
] 18(13 11( 5 10 13} 2 13) 2 5) 3 13} 
— 3 15 (21 38 (21 20 41) 3 $3} le 1(16 
.3- „U 29 (38 63 (48 32 (41 6 (37 Su ü {26 
e4=- „5 SS ler 56 (7? 23 (55) 3155) 6 (47 SO 
do .6 23 (66 35 (87 324 72) 7 (74 13 (82 10 (66 
TO= .7 28(82 15 ($3 16 (81 5 37) 2 (89) 7 (84 
7-7 .8 20 (94 9497 15 2) 3695 295) 2132) 
NI- .9 SOT Diy So 111% 53250 2 1921 11.9.5 
293=-7.0 E(100) 2¢100) 8ť100) 2(100) 0(190 2119») 


The Brier P-Scores for the three individual tests aná 


also th2 BPS for the sample climatology and for chance art 


(b 


enen in Teole XYII which indicates that each of the 
three methods is superior to climatology and chance a* 
all forecast tines. This is encouraging and is alse an 


Meat on that the SCF approach of describing the synoptic 


nom Cing is a valid approach that merits U aher 
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RABE ESA VIL 


P-Score Statistics for 36, 54 and 72 h Dependent and 
Independent Data Sets After "Filtering of the Predictand 
Score of 0.0 is Perfect) 
Dependent Set Independent Set 
PERS EOFS ESSE PERS EOFS P&E 
Bo h .308 ER] «39/0 2454 383 
climo 902 .u90 
chance -500 .500 
ieh «367 .409 . 340 . 401 393 . 368 
climo . 500 2504 
chance 2 509 .500 
man .389 153 278 a .441 s417 
pmo .500 . 500 
Chance 25010 215/90 
investigation. Neither PERS ər EOFS alone show as much 
skill as P&E at all fcrecast ‘times. This ic truüe Lor bet: 
the dependent and independent data sets. The difference 
between the two data samples is shown by the large 
Meqgzadesicon ci the BPS at 36 h for the independent sampie in 
all three tests. This may again indicate na the 
independent sample is tco smali. The scores in Table XVII 
indicate the decrease in skill as time increases. IR 
believed that this decrease in skiil with respect to tims 
could be lessened to some degree by the inclusion of 
ERES UNostic EOF islas. This was indicated by Shaffer (1982) 


as one of the 


= 


TOT 


largest drawbacks o 


Hh 


tre EOP tropical motion 


ecasting procedure. 





The third and simplest method of skill analysis is the 
percent correctly forecast (PCF) method. Table XVIII lists 
the ECFs ottained for each test for both the dependent and 
independent data sets. This table again illustrates the 


superiority of the P&E test after tha initial 36h time 


period. Notice again for the small independent data ser at 
54 h, the EOF method is superior to bot of the other 


methcds which is believed to be due to the clustering 
ERrect. 


TABLE XVIII 


pas 


Percent IN Forecast Events from tha Regression Metho 
Defining Recurvature as a Predictand Exceeding the 
Climatology Probabilities of Recurvature 
Score of 100 is Perfect} 


SD a’ A A A A A A A AO TT A A AM UND A A = =e ———À cA— 4D O LS A O A A ED EP au AA A A A A A OY A A A A A AY A A A A A A A A A A A A HP 


Dependent Set Independent Set 
PERS BORS PEE PERS EQUES PEE 
Boch ni. 7020 16% 3 Wile? OS Ze 
54 h nb. 3 7033 JN? p2.1 26.7 02.1 
Eh puso 623283 2s 7 73. 1 68.4 70.3 


This sinpie method of testing skill was included as an 


attempt to ccmpare the three tested procedures  adains- the 
p p Į d 


ip 


models currently being used by JTWC Guam for guidance. TE 
discussion and comments pertaining to the results of the 


model comparisons are included in Appendix A. 
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Elsberry (1979) discusses the inhsrent problems and 
advantages of dynamic tropical cyclone models. He states 
that the lack of vertical and ac one resolution in the 
models, as well as a poor understanding of The environmental 
processes that occur in the vicinity of a tropical storm, 
are limiting the progress toward more accurate forecasts. 
Although nested grid models (Harrison, 1981; Harriscn and 


haeor ino, 1982) have been usai to improve the resolution, 


(D 


there continue to be problems of data acquisition on th 


scales of interest. This problem has no scluticn in th 


(D 


ctr 
TF 
D 


immediate or extended time frames, which necessitates 
need for other reliable and less ccstly methods. 
The application cf a linear regression scheme to a 


Statistical model usirg a combination of persistence-related 


parameters and an EOF representaticn of th? synoptic forcing 
was completed. The Veembinatiscn ) of thease two types of 


parameters have shown to De usuper-or for forecasting 


tropical cyclone recurvature as defined. 
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Inem ca tub rDuUt ss | OL the EOF approach to forecasting 


tropical cyclone recurvature and motion (Shaffer, 1982) 


again be stressed: 


e The method is extremely inexpensive; 


(D 


must 


e The forecasts can be done ona handheld programmable 


GalCulatcr: 


e The method has been shown to be more accurate 


ESE 


forecasting recurvature than persistence and synoptic 


parameters alone; and 


e The method has shown to give superi 
cenparedi to mene Offic al. JTWC for 
analysis of the 500 mb fialds (Shaffer 


wo 


It is believed that the inclusion of prognostic synooti 


fields, as suggested by Shaffer (1982), would improve 


- fcrecasts when 
casts Using ens 
1982 


“m = P | là - eds - j ta o 
percentage of the explained variance of the predictand 


This would also be reflected aS giving batter  rzssul*ts 


the recurvazture study described here. 


An independent oferaticnal evaluation of the schemes 


presented here needs to be done to determine the validi 


the method. 


Meonder-n. £Eion of recurvature needs t5 be modified. 


tropical storm moving 314? is considered 


non-recurvature case where a storm moving 3169 is cons: 


a recurvature case. PEO = 2s OL E=scurvatuze ba 


the departure of tne tropical cyclone from a forecast 


(D 
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ty 


of 
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or a straight line could be utilized instead cf using a 
binary definition as done in this research. 

It is understood that the application of new and 
superior methods of forecasting tropical cyclone recurvature 
and motion should be questionei. It is believed chat this 
change of methodclogy (ie, using an EOF representation of 
the synoptic fields), if adopted, wiil inprove the image of 
the military weather services by being a positive pioneering 


effort in an area that has be=n stagnant for many years. 





APPENDIX A 


MODEL COMPARISONS : 


À. FURPOSE 


The purpose of this section is to supply the Western 


pa 


Peemric tropical cyclcne forecaster with information on the 
Meer ous Guidance techniques for forecasting “ropical cycions 
Mei and, 22) Pare i cular, thelr teliability in forecasting 
recurvature. The defini on of» recurvature (Leftwich, 
1978,1982) is similar as described earlier except for the 
northeastward motwen during the 24 h (vice t2) prior to 
the end of the forecast period. This nol f cation wes done 


to accomodate the forecast warning verification times which 


mee 24, 48 and 72 h. 


E UBVGtU-e iS a net displacement northward of 315° 
meng the forecast i$ period Or ene. atea Deng ere 
northeastward notion during the 24 hours prior +9 tae end 


of the forecast period. 


B. THE DATA 

The data analyzed were assembied by Dr. Ted Tsui of <ıe 
EE Environmental Prediction Research Facility  (NEPRF) 
during the 1979-1981 Northwest Pacific typhoon seasons. The 


Gata were the best track information from the tropica 





a 


disturbance stage thrcugh the entire lifs cycle of each of 
the 28 storms or 1979 and 1980 and the 29 storms of 1981. 
Also included were the 24, 48 and 72 h positions of 12 
tropical cyclone forecast aids available at JTWC. Seven of 
these aids and the Joint Typhoon Warning Center official 
forecasts are evalvated below with respect to the 


recurvature decision. 


C. MODEL DESCRIPTIONS COMPARISON TESTS 


cf 
Sc 
iD 


The f£ellowiag paragraphs are descriptions SE 


3 


E 
(D 


ob jective forecast aides currently being used at JTWC. 
descriptions have been taken from the 1981 Annual Tropical 
Cyclone Report. All JTWC objective techniques, except those 
involving dynamic modeis, ars executed using operational 
data available prior to warning times  (00GMT and  12GMT). 
These dynamically-based objective aids area usually recsived 
EMEUU D wtho;n 2 to 5 hours Of a Specific warning time. 

An analogue methcd, abbreviated TOPS, provides three 
movement and intensity forecasts. These three forecasts are 
BOLT analogues o* straight-moving cyclones, of recurving 
mees2ones and a combination or analogues including straight, 
Bo@urving, and all other cyclones that do not meet the 


Epecnfic cirteria of those t4o categories. 





An 2xtrapolation technique, abbreviated XTRA, provides a 
track connecting the 12-h Sila preliminary best track 
position and the current warning position, which is then 
extrapolated to 24 and 43 h. HPAC, which is a combination 
of XTRA and a climatological aid, provides a forecast of the 
track and intensity changes based on the initial position of 
the system. The ciimatological data base includes the 
period from 1945 to 1973. The (BPAC) forecasts are based on 
Anding the past motion of the tropical cyclone with the 
climatology forecast positions. This program requires cnly 
ae xas Instrument (T1-59) calculator +o Jsmerate 12 to 72 h 
forecast positions. The blending routine gives less weight 
to persistence at each succeeding forecast interval. 

An updated version of the HATRACK/MOHATT steering 
program (now called CYCLOFS and designated as cCY50 for the 
500 mb level) can provide steering forecasts at many 
different levels. The program can be run in the modified 
(includes a 12-h persistence bias) or unmodified versions 
applied to either analysis or prognostic fields. The 


program advects a point vortex on a pra-seiected analysis 


fps 
ho 
(D 
< 
iD 
| 
U) 
l i 
$, 
Or 
a 
Hm 


and/fer smoothed prognostic field at designaté 


time steps zhrough 72 h. In the modified version, che 
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program uses the previous 12-h history position to compute 
the 12-h forecast errcr and applies a bias correction to the 
forecast positions. In 1981, only the modified version, in 
the prognostic mode for the 500 mb level was verified. 

.The dynamic Tropical Cycione Model (TCMO) is a coarse 
mesh (220 km) primitive equation model. The digitized 
Sverone warning position is bogused into the 850 ab wind and 
temperature fields of the FNOC global band analysis. 
Hemispheric mccel forscast data are used on the boundaries. 
Two versions are currently run: One runs from forecast 
fields and is available on an automatic request basis, and 
@ae Orne r "zs initialized on the analysis fields notsd 
earlier. The model based on the analysis fields was tested 
mae this study. 

The NTCM or "nested" primitive equation tropical cyclones 
model is initialized on FNOC 12-h forecast fields. The 
Medel covers a limited, but relocatabplse, tropica. domain 
Meme three layers in tne vertical. Tha fines Scale, OT 
inested! grid, covers a 1200 km Square area with a 41 kn 
Ra. Spacing and moves to keep a 850 mb vortex in its 
center. This model is integrated with a coarse channel 


model grid with a grid spacing of 205 km over a 6490 oy 4700 


ON 
Ul! 





km demain. once initialized, the model runs independent of 
the remaining FNOC forecast fields. The NTCM is available 


via automatic request for 00 and 12 GUT forecast fields 


Only. 


Dee SCREENING PROCEDURE FCR MODEL TESTING 


Several differert requirements were taken nO 
Sorstderaticon Comm inclusion in this comparison study that 
were not considered fcr the ZOP analysis scheme. Tnclusson 


was based on: 


° eag the ‘best track! past zonal and meridional 
G@esprde=stents (after the initial 12 h of the track) 
vice the past warning dispiacements usec in the EOF 
study to test for eastward displacement (mer 
considered case) ; 

e Pseudo-irdecendence is again assumed by requiring time 
separaticns of 36 à within the same storm; 


e A 30 h past history was not required; 


ee meeroee cal Cyclones wera Sencluded,  zather than just 
promos storms and -vphoons; ard 


EE. 0:55 nor-n or^ ene equator, vise 109 north, were 
PG rude G « 
Table XIX lists the 24, 48 and 72 k r ecurvaturs 


probabilities for the three year sample. When compared to 
the corresponding values of the sample climatolcgy listed in 
Table II, the difference between “he inclusion criteria 


become evident. The inclusion of the tropical depressions 
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in the area south of 109 north is tne largest contributor to 
the lower values of recurvaturs probebility. A high annual 
variability is indicated in this small sanpie. There is 
generally a slow increase of probapility of Tecurvature with 


time as was the case in Table II. 


TABLE XIX 
EIA E ETS SEE 
mc A mu en.) een eee 
1979 1980 1981 

tau rcv nrcv prob FEV” DEeevy prob ECY nicy prob 
24 39 94 298 36 78 2.3105 26 73 «203 
48 41 84 2379 35 55 350 23 63 2:297 
72 25 81 226 29 DES . 354 21 52 Ze 


The HSS for the three separate years ır= listed in Table 
RX. The HSS in this data set remains as described earlier 
ween a slight modification. This difference is that the 
n fication of che model's predicted recurvature is pasad 
entirely upon the predictands matching the actual 
recurvature criteria as defined. There is no climatological 
threshhold that is set for inclusion as a recurvature or 


non-recurvature as was the case for the 


(D 


arlier data seupie. 


Again, the most noticable feature of the tabie is the 


(1 


general decrease of skill with respect to tine. The 
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exception to this generality is the BPAC model performance 
for 1979 and the TCMO model for 1981 where the HSS increased 
with time. An interesting comparison can be shown when 
comparing JTWC with the TCMO forecast aide for 1981. There 
is a definite correlation betwean TCMO ani JTWC. This is an 
w cation of JTNC Ecliowing the guidances of this dynamic 
model quite closely in the extended time frames. 


Table XXI lists the BPS for the thr 


qi 


e time periods of 
interest. The BPS forecast predictands in this case are not 
on the interval between 1.0 and 0.0, but are identically 
caual 921.0 or 0.0. Any knowledge ož a comparison of BPS 
results using a continuous. predictand as compared to a 


binary predictand schene is available. Comparison ot wıe 


I 


Connon tame period of 72 betwesn the two entirely 
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independent data sets cf dif 
little meaning. Themers Of Gen TCMO model in 1981 sed BPAC 


meme 1979 are the only models that approach “he BPS for ths 


combination o£ SOFS and Persistence in the earlier study. A 
general trend toward lazger BPS are prevalent for ali years 
for most of the forecast aids. 

ompa ison of the PCF valuas indicated in Table XXII, 
and the HSS reveals the similarity of tha two skill scores. 
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Using the skili scores described in Chapter V for the 
1979-81 typhoon season, ncdel performance indicates a large 


Aia varaation of skill between the models. 
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